Measurements of if)(2S) decays to octet baryon-antibaryon pairs 
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With a sample of 14 x 10 6 ip{2S) events collected by the BESII detector at the Beijing Electron Positron 
Collider (BEPC), the decay channels ip(2S) — > BgBs (pp, AA, E°E°, E~H + ) are measured, and their branching 
ratios are determined to be (3.36 ± 0.09 ± 0.25) x 10 -4 , (3.39 ± 0.20 ± 0.32) x 10~ 4 , (2.35 ± 0.36 ± 0.32) x 10~ 4 , 



1. Introduction 

The branching ratios of ip(2S) decays into octet 
baryon-antibaryon pairs were measured by the 
BES-I and CLEOc collaborations, and the results 
differ significantly, as shown in Table 1. It is 
therefore important to make new measurements 
to help clarify these differences using the sample 
of 14x 10 6 xp(2S) events collected by BESII, which 
is the world's largest e + e~ ip(2S) sample. 



are several measurements for a for J/ip — > pp, 
and the recent result of a = 0.676 ± 0.036 ± 0.042 
given by BES [6] is quite close to Carimalo's pre- 
diction a = 0.69 [5]. However, there is only one 
measurement for tp(2S) — * pp, made by E835 [7]. 
Results for tp(2S) — ► pp are summarized in Ta- 
ble 2. The i/j(2S) -> pp events in BESII allow 
the measurement of a, which can be compared 
with the existing result and used to test hadron 
helicity conservation. 



Table 1 

Branching ratios of ip{2S) -> B 8 I? 8 (x 10~ 4 ). 



Channel BES-I [1] 



CLEO-c [2] 



pp 2.16 ±0.15 ±0.36 2.87 ±0.12 ±0.15 

AA 1.81 ±0.20 ± 0.27 3.28 ± 0.23 ± 0.25 

1.2 ±0.4 ± 0.4 2.63 ±0.35 ±0.21 

0.94 ±0.27 ± 0.15 2.38 ± 0.30 ± 0.21 



According to the hadron helicity conservation, 
the angular distribution of ip(2S) — » Bg,Bg, can be 
expressed as: 

dN 



oc 1 + a cos 



(1) 



dcos9 

where 9 is the angle between Bg, and the beam 
direction of the positron in the center-of-mass 
(CM) system. In the limit of infinitely heavy 
charm mass, hadron helicity conservation implies 
a = 1 [3] for both J/ip and tp(2S) decays to octet 
baryon anti-baryon pairs. 

Values of a for J/tp, ?p(2S) — > pp have been 
predicted theoretically based on first order QCD. 
In the prediction of Claudson, Glashow, and 
Wise [4] , the mass of the final baryon is taken into 
account as a whole, while the constituent quarks 
inside the baryon are taken as massless when com- 
puting the decay amplitude. In the prediction by 
Carimalo [5], mass effects at the quark level are 
taken into consideration. Experimentally there 



Table 2 

Predicted and measured values of a for ip(2S) 
pp. 



a value 


Source 




Predicted value: 


a = 0.58 


Claudson et al. 


[4] 


a = 0.80 


Carimalo [5] 




Measured value: 


a = 0.67 ±0.15 ±0.04 


M. Ambrogiani 


et al. [7] 



BESII is a large solid-angle magnetic spec- 
trometer which is described in detail in Ref. [8]. 
The momentum of charged particles is deter- 
mined by a forty-layer cylindrical main drift 
chamber (MDC) which has a resolution of 
o-p/p=1.78%-v/l +p 2 (p in GeV/c). Particle iden- 
tification is accomplished using specific ioniza- 
tion (dE/dx) measurements in the drift cham- 
ber and time-of-flight (TOF) information in a 
barrel-like array of forty-eight scintillation coun- 
ters. The dE/dx resolution is a dE/dx — 8.0%; 
the TOF resolution for Bhabha events is oto f = 
180 ps. Radially outside of the time-of-flight 
counters is a 12-radiation-length barrel shower 
counter (BSC) comprised of gas tubes interleaved 
with lead sheets. The BSC measures the en- 
ergy and direction of photons with resolutions of 
a E /E ~ 21% /VE (E in GeV), o> = 7.9 mrad, 



2 



3 



and a z — 2.3 cm. The iron flux return of the 
magnet is instrumented with three double layers 
of proportional counters that are used to identify 
muons. 

Monte Carlo (MC) simulation is used for mass 
resolution and detection efficiency determination. 
In this analysis, a GEANT3 [9] based MC package 
(SIMBES) with detailed consideration of the de- 
tector performance (such as dead electronic chan- 
nels) is used. The consistency between data and 
MC has been carefully checked in many high pu- 
rity physics channels, and the agreement is rea- 
sonable [10]. 

The data samples used for this analysis con- 
sist of 14.0 x 10 6 (1 ± 4%) ^(2S) events [11] and 
6.42(1 ± 4%) pb _1 of continuum data at s/s = 
3.65 GcV [12]. The decay channels investigated 
are ip{2S) -> pp, AA, and E~E + , where A 

decays to ^"(63.9%), E° decays to A 7 (100%), 
and S" decays to A7r~(99.9%). 

2. Event selection and branching ratio de- 
termination 

2.1. ip(2S) pp 

The experimental signature for the decay 
-0(2S') — ► pp is two back-to-back, oppositely- 
charged tracks, each with a momentum of 
1.586 GeV/c. The main backgrounds are: 
Bhabha and dimuon (e + e~ — > events, 
H2S) -► e+e~, n+fi-, tt+tt", K+K~, ip{2S) -► 
7Xcj(j=o,i,2) -> jtt + tt-(K + K~, pp), ip(2S) -> 
ir°ir°J/ij -> 7r 7r°e+e-( M + ^-), ^(2S) -> 7T°p|>, 
etc. 

The event selection requires two well recon- 
structed and oppositely charged tracks. Each 
track is required to be well fitted to a three 
dimensional helix, be in the polar angle region 
| cos#| < 0.8, and have a momentum greater than 
70 MeV/c in the xy-plane. The point of clos- 
est approach of each track to the bcamline is re- 
quired to be within the interaction region which 



is defined to be ±20 cm longitudinally and 2 cm 
radially. 

In order to remove cosmic rays, the difference 
between the time-of-fiights of the positive and 
negative tracks, \t + — is required to be less 
than 4.0 ns. Protons and antiprotons are required 
to be identified by the TOF; the measured time- 
of-flight of the track must be closest to the predic- 
tion for the proton/anti-proton hypothesis. Since 
tp(2S) — » pp is a back-to-back two-body decay, we 
require the acollincarity angle of two tracks to be 
less than 5°. The deposited energy in the BSC 
of the positive particle is required to be less than 
0.75 GeV to remove possible e + e~ final state con- 
tamination. Finally, the energy sum (calculated 
from the track's momentum) of the two tracks is 
required to be within 130 MeV of the expected 
sum, 3.686 GeV, and the momentum of the neg- 
ative track is required to be within 150 McV/c of 
the expected momentum 1.586 GeV/c. 

Events surviving the selection criteria are 
shown in Fig. 1 as dots with error bars. The 
same selection criteria have been applied to back- 
ground events generated by the MC and nor- 
malized according to branching ratios listed in 
PDG(2006) [13], and 38.1 background events sur- 
vive and are shown as the dashed line in Fig. 1. 
The data are fitted by a MC histogram for the sig- 
nal plus a background function which corresponds 
to the 38.1 simulated background events and a 
flat distribution to describe the remaining back- 
ground. From the fit, the number of pp events is 
determined to be 1618.2 ± 43.4, where the error 
is statistical. 

2.1.1. Angular distribution of pp 

To obtain the parameter a for iJ}(2S) — ► pp, the 
cos 8 dependence of the event selection efficiency 
must be taken into account, which is determined 
using a flat angular distribution (a=0) in the MC 
simulation; see Fig. 2(a). However, there are im- 
perfections in the MC simulation, which will dis- 



4 




1 .440 1 .540 1 .640 

Pproton (GeV/c) 



Figure 1. The fitted proton momentum spectrum. 
The dots with error bars are data, the histogram 
is the fit to the data including the signal shape 
from MC and all backgrounds, and the dashed 
line is the background. 



tort the efficiencies determined by the MC as a 
function of cos 9. In order to reduce this system- 
atic error, a correction to the MC efficiency is 
made [6]. The correction factor f c (cos9) is de- 
fined as: 

^Data , 



/ C (COS0) 



SMC 



3S0)=n ££ ^(cOS0), 

\ L e M c(i) 



where i denotes the selection criterion, £Data(i) 
is the efficiency determined for data for criterion 
i, and EMc{i) is the efficiency from the MC for 
criterion i. The corrected MC efficiency is then: 

e' mc (cqs0) = £ MC (cos 6) x / c (cos0). 

Due to the limitation on the number of 
-0(25*) — > pp events, the "reference" channel 
J ftp — > pp is chosen to determine the correc- 
tion factor due to its higher statistics and sim- 
ilar kinematics. The selection criteria related to 



the energy and momentum for ip(2S) — * pp are 
scaled to the reference channel J/tp — > pp. Then 
following the re-weighting procedures in Ref. [6] 
for our selection criteria, the correction function 
f c (cos9) is obtained and is shown in Fig. 2(c). 
With emc (cos 8) denoting the efficiency obtained 
from ip(2S) — ► pp MC and f c (cos 9) the correction 
function for the efficiency, we fit the measured an- 
gular distribution of ip(2S) — > pp data with the 
function N(cos9), 

N(cos9) = 7V x(l+acos 2 6) x smg (cos 6) x / c (cos 6), 

as shown in Fig. 2(d). The fit uses a binned 
X 2 minimization method in the angular range 
cos (9 G [-0.7,0.7] and gives xLn = 10 - 88 for 12 
degrees of freedom. The fitted value of the pa- 
rameter a is 0.85 ± 0.24, where the error is sta- 
tistical. 

As a consistency check, we also obtained 
f c (cos9) directly from the ip(2S) — > pp sample, 
and the fitted result obtained using this correc- 
tion yields a = 0.83 ± 0.24, but its systematic 
uncertainty is 0.14, mainly due to the lower statis- 
tics of the ip(2S) — ► pp sample, and much larger 
than the systematic error of 0.04 determined us- 
ing f c (cos9) obtained from the J/iJj — > pp sam- 
ple (see section 3.1). This demonstrates that 
f c (cos9) determined from the J/ip — > pp sam- 
ple improves the systematic error on a without 
changing its central value and statistical error. 

2.1.2. Branching ratio of xp(2S) — ► pp 

The selection efficiency is determined using 
1 x 10 5 V(2S) -► pp MC events. The MC- 
determined efficiency is emc = (34.4±0.2)%, and 
the branching ratio is determined to be: 

Br{ip(2S) -> pp) = (3.36 ± 0.09) x 10~ 4 , 

where the error is statistical. 

2.2. i/)(2S) -> AA 

Candidate events require four well recon- 
structed charged tracks. The positive (negative) 
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Figure 2. (a) The selection efficiency versus cos 9 
obtained from MC; (b) angular distribution of 
background events from MC, which survive the 
same selection criteria as used for data; (c) the 
correction obtained from data (/ c (cos0)) to the 
MC efficiency; and (d) the angular distribution 
of candidate tp{2S) — > pp events. 

charged track with the higher momentum is as- 
sumed to be the proton (antiproton); the other 
two are assumed to be the n + and n~ . The two 
pir pairs are required to pass the A's vertex find- 
ing algorithm successfully, and the sum of the A 
and A decay lengths must be greater than 0.02 
m (see Fig. 3). The sum of the A and A ener- 
gies must be in the region from 3.60 GeV to 3.81 
GeV (sec Fig. 4). The missing momentum of the 
events should be less than 0.18 GeV/c. and the 
difference between the measured mass of Mp^+ 
and its expected value. Ma, should be less than 
12 Me/c 2 (three times the resolution of the Ma). 

The events that satisfy the selection criteria are 
shown in Fig. 5 as dots with error bars; they are 
fitted by a histogram of the signal shape from 



Figure 3. The sum of the A and A decay lengths. 
The histogram is the signal shape from the MC 
plus simulated background, the dots with error 
bars are data, and the shaded histogram is the 
background. 



MC plus a background function which describes 
the simulated backgrounds and a flat distribution 
to describe any remaining sources. The simulated 
backgrounds are mainly from ip(2S) — > S°I] and 
ip(2S) — ► AS + c.c. and normalizing according to 
branching ratios from PDG(2006), a total of 32 
background events arc obtained. The final num- 
ber of signal events from the fit is 337.2 ± 19.9. 

The ip(2S) — > AA — > ppir + -K~ efficiency is de- 
termined to be e MC = (17.4±0.2)% using 2 x 10 5 
MC-simulated signal events. The branching ratio 
is then: 

Br{^{2S) -> AA) = (3.39 ± 0.20) x 10~ 4 , 
where the error is statistical. 
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Figure 4. The sum of the A and A energies. The 
histogram is the signal shape from the MC plus 
backgrounds. The dots with error bars are data, 
and the dashed line is the main background from 
ip(2S) -> The peaks at 3.1 and 4.2 GcV 

are from ip(2S) -> 7r 7r°J/V>, J/ip -> AA and 
ip(2S) -> Tr+n-J/iP, J/i< -> e+e- (or m+m"), 
respectively. 

2.3. -0(25') -»■ 

Candidate events are required to have four 
well reconstructed charged tracks plus at least 
two good photons. The A and A are selected 
using the method described in Section B. The 
missing momentum of the events is required to 
be less than 0.25 GeV/c. The \ 2 of the four 
constraint (4C) kinematic fit to the hypothesis 
ip(2S) — > pp7r + 7r - 77 must be less than 20. The 
difference between the measured mass of Mp 7T + J 
and its expected value, M^o , should be less than 
36 MeV/c 2 (three times the Mgo resolution). 

The events that survive selection are shown 
in Fig. 6 as dots with error bars; they are fit- 
ted by a histogram of the signal shape from MC 



Figure 5. The fitted A mass spectrum. The dots 
with error bars are data, the histogram is the fit 
to data which includes the signal shape from MC 
plus all backgrounds, and the dashed line is the 
background. 



plus a background function which describes the 
simulated backgrounds and a flat distribution to 
describe the remaining background. The main 
backgrounds are from ip{2S) — * AA, ip(2S) — > 
7X_cj(./=o,i,2) -> 7AA, iP(2S) - S°H°, ^(2S) - 
AE°+c.c. and i/)(2S) -> S°S°+c.c., and normaliz- 
ing using branching ratios from PDG(2006), 16.5 
background events are obtained. The final num- 
ber of signal events from the fit is 59.1 ± 9.1. 

The ip(2S) -> AA77 -> ^+71-77 

efficiency is determined to be tuc = (4.4 ±0.1)% 
using 2 x 10 5 MC generated signal events. The 
branching ratio of signal channel is then: 

Br{i)(2S) -> S°S°) = (2.35 ±0.36) x 10~ 4 , 

where the error is statistical. 
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Figure 6. The fitted S° mass spectrum. The dots 
with error bars are data, the histogram is the fit 
to data which includes the signal shape from the 
MC and all backgrounds, and the dashed line is 
the background. 



Figure 7. The S - and S + energy sum. The his- 
togram is the signal shape from MC plus back- 
ground. The dots with error bars are data, and 
the shaded area is the sum of simulated back- 
grounds. 



2.4. -0(25) -> H-H+ 

Candidate events require six well reconstructed 
charged tracks. The positive (negative) charged 
track with highest momentum is assumed to 
be the proton (antiproton); the other four are 
assumed to be 7rs. Looping over all possi- 
ble pn~ , pTt + combinations in an event, the 
one which successfully passes the vertex find- 
ing algorithm and has the smallest value of 
■s/iMpn- - M A ) 2 + (M pw+ - M A ) 2 is selected for 
further analysis. The energy sum of the S~ and 
S + should be between 3.593 and 3.779 GeV (see 
Fig. 7) , and the missing momentum of the events 
should be less than 0.15 GeV/c. The difference 
between the measured mass of M™+ w + and its 
expected value, M^+, should be less than 18 
McV/c 2 (three times the Mb- resolution). 

The events surviving selection are shown in 



Fig. 8 as dots with error bars, and they are fit- 
ted by a histogram of the signal shape from MC 
plus a background function which describes the 
simulated backgrounds and a flat distribution to 
describe remaining background. The main back- 
ground is from ip(2S) — ► Tt + n~ J /ijj — > 7r + 7r _ AA, 
and normalizing by PDG(2006) branching frac- 
tions, 11.3 background events are obtained. The 
final number of signal events from the fit is 
67.4 ±8.9. 

The V(2S) S"S+ -> AAtt+tt- -> 

pp7T + 7r~7r + 7r~ efficiency is determined to be 
£mc = (3.9 ±0.1)% using 2 x 10 5 signal events 
generated by MC. The branching ratio of the sig- 
nal channel is then: 

Br(^(2S) -> E-E+) = (3.03 ± 0.40) x 10~ 4 , 

where the error is statistical. 
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Figure 8. The fitted S~ mass spectrum. Dots 
with error bars are data, the histogram is the fit 
to data which includes the signal shape from the 
MC and all backgrounds, and the dashed line is 
the sum of the backgrounds. 



3. Systematic error 

3.1. ■0(25') — > pp angular distribution 

The systematic error on a in i^{2S) — > pp decay 
from the tracking reconstruction is determined 
using different MDC wire resolution models in 
the MC simulations, which changes a by 2.7%. 
When the fit parameter of the efficiency correc- 
tion curve / c (cos#) is changed by la, a changes 
by 2.3%. The performance of the BES detector 
has small differences between the time when the 
58 x f 6 J/ip events were obtained and when the 
14 x 10 6 ip(2S) events were obtained. Using pa- 
rameter files describing the performance of BES 
detector at these two data taking periods, the ef- 
fect of this variation on a is determined to be 
2.2%. The effect of the background uncertainty 
on a is negligible. Adding these contributions in 



quadrature gives a total systematic error of 4.2%. 

3.2. Branching ratios 

The systematic errors on the branching ratios 
are mainly from the uncertainties in the MDC 
tracking, a, the hadronic interaction model, back- 
ground estimations, and differences between data 
and MC for the A vertex finding, decay length 
requirement, and kinematic fitting. 

The MDC tracking gives a systematic error of 
about 2% for a proton or anti-proton [10] and 1% 
for a low-energy tt, which is determined from the 
channel ip(2S) — ► tt + tt^J/i/j — > it + it~ p + pT . The 
detection efficiency depends on the angular dis- 
tribution of the baryon pair. For pp decay, when 
changing the a value by lcr, the branching ra- 
tio changes by 2.4%; in the other three channels, 
a = 0.5 is used as a nominal value, the maximum 
differences for |e Q=0 .5 - £q=o| and |e Q=0 .5 - e a =i| 
are taken as systematic errors, they are 6.5%, 
7.6%, 6.8%, respectively. The uncertainties of the 
detection efficiencies caused by assumed flat an- 
gular distributions for secondary decay of baryons 
are much smaller than those from angular distri- 
butions of ip(2S) to baryon pair primary decays, 
and are therefore neglected here [14]. Different 
simulation models for the hadronic interaction 
(GCALOR/Geant-FLUKA) [15,16] give different 
efficiencies, giving systematic errors of 2.18%, 
0.46%, 0.00%, 1.08% for the studied channels, re- 
spectively. The background uncertainty is stud- 
ied by changing the nominal branching ratios of 
the backgrounds which have large statistical er- 
rors. If the branching ratios of the background 
channels are changed by 100% in the pp, AA, 
and channels, the changes in the branch- 

ing ratios in the signal channels arc 0.1%, 1.0% 
and 0.2%, respectively. For the E°E° channel, 
where the shape of the simulated backgrounds is 
in good agreement with the data in the A7 in- 
variant mass distribution, the branching ratios of 
backgrounds are only changed by 20%, resulting 
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in a change of the branching ratios of the signal 
channel of 2.3%. According to the reference chan- 
nel J/ip — > AA [17], the secondary vertex finding 
of A gives a systematic error of 0.7% for each A 
vertex, and the requirement on the sum of the 
decay length contributes 1.4%. 

In the branching ratio determination of four 
channels, the continuum contribution must be 
subtracted. The continuum data are also selected 
with the same criteria as for the ijj(2S) decay 
signal channels, and the number of the surviv- 
ing events times a luminosity normalization fac- 
tor is taken as a systematic error. The kinematic 
fit of pp7r + 7r~77 in ip(2S) —* gives a sys- 

tematic error of 7.6% from the reference channel 

1p(2S) -> TT+TT-J/ij -> Tr+TT-TT+TT-TT [18]. The 

uncertainty on the total number of ip(2S) events 
is 4%. The systematic errors of the acollinearity 
angle, E B § region, baryon mass (or momentum), 
and Pmiss requirements are studied with corre- 
sponding J/tp — > pp decays. 

In the ip(2S) — > pp selection, the systematic er- 
rors due to the uncertainties from particle iden- 
tification, the cosmic ray veto, and the deposited 
energy criterion are studied by this channel itself. 
All the systematic errors in the branching ratio 
measurements are summarized in Table 3. 



ip(2S) -► pp is measured to be 0.85 ± 0.24 ± 0.04, 
which is in agreement within la with the E835 
result [7], and close to Carimalo's prediction [5]. 

Table 3 

Systematic errors in the branching ratio measure- 
ments (%). 



Source 


pp 


AA 






MDC tracking 


4 


4.5 


4.5 


5.7 


PID 


2.4 








Cosmic Ray Exc. 


0.9 








Deposit Energy 


0.9 








Acol. angle 


0.9 








Vtx. finding 




1.4 


1.4 


1.4 


Decay length 




1.0 


1.0 




E BBl M B (or P B ) 


0.8 


0.6 


1.6 


1.6 


1 rniss 




1.6 


0.5 


1.7 


7 tracking- 






4 




Kinematic fit 






7.6 




Bg. Esti. 




1.0 


2.3 


0.2 


Continuum data 


0.8 


1.0 






a value 


2.4 


6.5 


7.6 


6.8 


Hadronic Interaction 


2.2 


0.5 




1.1 


N 4>(2S) 


4 


4 


4 


4 


Total error 


7.3 


9.4 


13.4 


10.3 



4. Summary and discussion 

Based on 14 x 10 6 ip(2S) events, the branch- 
ing ratios of ip(2S) -> pp, AA, E S°, and 
are measured, the results are listed in Table 4, 
together with the ratios of ip(2S) — > BB to 
J ftp — > BB. They are in agreement with the 
results published by the CLEO collaboration [2] 
within 2a for pp and within la for the other three 
channels. The differences of the branching frac- 
tions between current measurements and those of 
BESI are 2.5c, 3.1cr, 1.5c, 3.5er for the four chan- 
nels, respectively. 

The angular distribution parameter a for 
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Table 4 

Branching ratios of ip(2S) decays into baryon 
anti-baryon pairs. The first error is statisti- 
cal and the second systematic. The value Q is 
BR(iP{2S) -> BB)/BR(J/tp -> BB). The J/ip 
branching ratios are taken from Ref. [6] for pp, 
Rcf. [17] for AA and and Ref. [13] for 



modes 


BRs (xlCT 4 ) 


Q (%) 


PP 


3.36 ±0.09 ±0.25 


14.9 ± 1.4 


AA 


3.39 ±0.20 ±0.32 


16.7 ± 2.1 




2.35 ±0.36 ±0.32 


16.8 ±3.6 




3.03 ±0.40 ±0.32 


16.8 ±4.7 
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